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We report the nature of transitions from nonsynchronous to complete synchronization (CS) state 
in arrays of time-delay systems, where the systems are coupled with instantaneous diffusive coupling. 
We demonstrate that the transition to CS occurs distinctly for different coupling configurations. In 
particular, for unidirectional coupling, locally (microscopically) synchronization transition occurs in 
a very narrow range of coupling strength but for a global one (macroscopically) it occurs sequentially 
in a broad range of coupling strength preceded by an intermittent synchronization. On the other 
hand, in the case of mutual coupling a very large value of coupling strength is required for local 
synchronization and, consequently, all the local subsystems synchronize immediately for the same 
value of the coupling strength and hence globally synchronization also occurs in a narrow range of 
the coupling strength. In the transition regime, we observe a new type of synchronization transition 
where long intervals of high quality synchronization which are interrupted at irregular times by 
intermittent chaotic bursts simultaneously in all the systems, which we designate as global inter- 
mittent synchronization (CIS). We also relate our synchronization transition results to the above 
specific types using unstable periodic orbit theory. The above studies are carried out in a well known 
piecewise linear time-delay system. 

PACS numbers: 05.45.Xt,05.45.Pq 



I. INTRODUCTION 

Numerical and experimental investigations of chaotic 
synchronization in coupled nonlinear systems have been 
receiving much attention in recent years. This phenom- 
ena is omnipresent and plays an iinportant role in diverse 
areas of science and technology In the synchroniza- 

tion process, two identical chaotic systems do not always 
necessarily synchronize perfectly. Rather, long intervals 
of high-quality synchronization are interrupted at irregu- 
lar times by intermittent chaotic bursts and such chaotic 
bursts along with the synchronization are called on-off 
intermittency [3] . It has been shown that on-ofF intermit- 
tency is a frequently occurring instability preceding typi- 
cal synchronization transitions in diverse dynamical sys- 
tems, mediated by unstable periodic orbits (UPOs) [J]. 
Further, as the coupling parameter is increased, a pe- 
riodic orbit embedded in the attractor in the invariant 
synchronization manifold can become unstable for per- 
turbations (such as noise and/or parameter mismatches) 
transverse to the manifold. This is called a bubbling bi- 
furcation, which leads to the formation of riddled basins 
of attraction in the invariant manifold inducing intermit- 
tent bursting (see Ref. [s^ for more details). There exists 
another type of bifurcation, called blowout bifurcation, 
induced by changes in the transverse stability of an in- 
finite number of UPOs. Among these UPOs, some are 
transversely stable and others are transversely unstable 
near the bifurcation. 

It is a well accepted fact that on-off intermittency is a 



common phenomenon which occurs in a wide variety of 
natural systems, including neural networks [EB], biolog- 
ical systems [1] , laser systems [1, [l^l , electronic circuits 



11 , |12| | , complex networks 11311 . coupled chaotic systems 
14 [, earthquake occurrence jl5j and other physical sys- 
tems such as Hamiltonian systems and self-driven par- 
ticle systems [l^, [l3l . Specifically, it has been reported 
that the dynamics of clusters in a network can exhibit 
an extreme form of intermittency [l8| : A substantial 
percentage of synchronized nodes forms a giant cluster 
most of the time, while many small clusters can also 
occur at other times. Thus the cluster sizes can vary 
in a highly intermittent fashion as a function of time. 
Recently, it has been shown [l^ that the transition to 
intermittent chaotic synchronization (in the case of com- 
plete synchronization (CS)) for phase-coherent attractors 
(Rossler attractors) occurs immediately as soon as the 
coupling parameter is increased from zero and for non- 
phase-coherent attractors (Lorenz attractors) the transi- 
tion occurs slowly in the sense that it occurs only when 
the coupling is sufficiently strong known as delayed tran- 
sition. 

It has been already shown that the transition from 
nonsynchronization to any type of synchronization is pre- 
ceded by intermittent synchronization in coupled chaotic 
systems. For example, intermittent lag synchroniza- 
tion (ILS) [13 J intermittent phase synchronization (IPS) 
(21I and intermittent generalized synchronization 
(IGS) [2 31 are some of the synchronization transitions 
characterized by the intermittent behavior as a function 
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of a coupling parameter. Recently, IGS has been numeri- 
cally observed in unidirectionally coupled time-delay sys- 
tems [13] ■ It has been found that the onset of generalized 
synchronization is preceded by on-off intermittency and 
the transition behavior is different for different coupling 
schemes. In particular, the intermittent transition occurs 
in a broad range of coupling strength for error feedback 
coupling configuration and in a narrow range of coupling 
strength for direct feedback coupling configuration be- 
yond certain threshold values of the coupling strength. 

The transition between various types of synchroniza- 
tion and their mechanism are not yet well understood 
especially in time-delay systems. Further, the dynam- 
ics of a large ensemble of coupled time-delay systems 
such as regular and complex networks are not yet well 
studied and only a very few studies are available in the 
literature (25l , [26| . The study of synchronization in en- 
sembles of time-delay systems has been receiving central 
importance recently in view of the infinite dimensional 
nature and feasibility of experimental realization of time- 
delay systems. Particularly, considerable attention is be- 
ing paid to time-delay systems with instantaneous cou- 
pling due to their extensive applications in different fields 
such as signal and image processing, pattern recognition, 
chaotic neural networks, secure communication and cryp- 
tography [23-113 ■ particular. In Refs. [I^IIO], it has 
been demonstrated that in chaotic communication exper- 
iments, the time-delayed optical fibre ring laser system 
is capable of transmitting the encoded signals with the 
speed of IGb/sec data range over a long distance fibre- 
optic channel (w 120Km). 

Motivated by the above, we will investigate here syn- 
chronization transitions in an array of coupled time-delay 
systems with different (instantaneous) coupling configu- 
rations. Particularly, in this paper, we demonstrate that 
the transition to CS occurs distinctly for different cou- 
pling configurations in a regular array of coupled time- 
delay systems. In an unidirectional array the transition 
from nonsynchronization to CS occurs locally (micro- 
scopically) in a narrow range of coupling strength and 
globally (macroscopically) the systems synchronize one 
by one with the drive system as a function of the coupling 
strength, which is known as sequential synchronization. 
But in a mutually coupled array, every individual sys- 
tem synchronizes immediately in a narrow range (after 
a large threshold value) of the coupling strength and so 
globally the synchronization transition is immediate as a 
function of the coupling strength in contrast to sequential 
synchronization. It is also to be noted that in the transi- 
tion regime we observe a new type of synchronization be- 
haviour called global intermittent synchronization ( GIS) 
where long intervals of high quality synchronization are 
interrupted by large desynchronized chaotic bursts simul- 
taneously in all the systems in the array. 

To understand the two distinct transition scenarios, 
we focus on the theory of unstable periodic orbits, which 
are the basic building blocks of chaotic and hyperchaotic 
attractors. The sequential and immediate synchroniza- 
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FIG. 1: (Color online) The probability of synchronization 
^i{e) as a function of the coupling strength e. The system in- 
dex i illustrates the occurrence of sequential synchronization 
transition to CS in unidirectionally coupled piecewise linear 
time-delay systems (Eq. HJ). The black color indicates the 
absence of synchronization ($i(e)=0.0), whereas the yellow 
color (light gray) represents the occurrence of CS ("I'i(£) = 1.0). 

tion transitions to CS are characterized by calculating 
the probability of synchronization and the average prob- 
ability of synchronization as a function of the coupling 
strength. The existence of intermittent synchronization 
is corroborated by using a spatiotemporal difference and 
a power law behavior of the laminar phase distributions. 

The remaining paper is organized as follows: In Sec. 
ITTl we will explain the occurrence of sequential synchro- 
nization preceded by intermittent synchronization in an 
array of unidirectionally coupled piecewise linear time- 
delay systems and in Sec. Illli we consider a mutual cou- 
pling configuration and explain the occurrence of instan- 
taneous synchronization transition in the array. Further 
we demonstrate the existence of GIS and provide a pos- 
sible mechanism for the occurrence of this new synchro- 
nization transition to CS in the array. Finally, we discuss 
our results and conclusion in Sec. IIVI 

II. SYNCHRONIZATION IN A PIECEWISE 
LINEAR TIME-DELAY SYSTEMS: LINEAR 
ARRAY WITH UNIDIRECTIONAL COUPLING 

We consider the following unidirectionally coupled 
time-delay systems of the form 

±1 = -axi{t) + I3f{xi{t - t)), (la) 
it = -axi{t) + I3f{xi{t - t)) + e[xi^i{t) ~ Xi{t)lllh) 

where i — 2,3, ■ ■ ■ ,N. We choose an open end boundary 
condition, a, /3 are system parameters, r is the time- 
delay and e is the strength of the coupling between the 
systems. The nonlinear function f{x) is chosen to be a 
piecewise linear function with a threshold nonlinearity, 
which has been studied recently [33| . 

f{x) = AF* - Bx. (2) 
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FIG. 2: (Color online) Snap shots of node vs node plots in- 
dicating sequential synchronization in unidirectionally cou- 
pled piecewise linear systems for different values of coupling 
strength, (a) e = 0.4, (b) e = 0.7, (c) e = 0.87 and (d) 
e = 1.1. 




FIG. 3: (Color online) The probability of synchronization 
^i{e) of selected systems (i = 15,20,30) and the average 
probability of synchronization (^(e)) in an unidirectionally 
coupled array (Eq. ((T])) as a function of the coupling strength 
e. 
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The system parameters for the piecewise linear system 
are fixed as follows: a = 1.0, /3 = 1.2, t = 6.0, 
A — 5.2, B = 3.5 and x* is the threshold value fixed 
at X* — 0.7. Note that for this set of parameter values a 
single uncoupled system exhibits a hyperchaotic attractor 
with three positive Lyapunov exponents (LEs) (see Ref. 
0). 

To demonstrate the nature of the dynamical transition 
to a complete synchronization regime, we consider an ar- 
ray of iV = 30 unidirectionally coupled identical piece- 
wise linear time-delay systems (each system hav- 



ing different initial conditions). Here, Xi{t) acts as the 
drive and the remaining systems {xi(t), i — 2,3, ■ ■ ■ , 30) 
as the response systems. In the absence of coupling 
[e = 0.0 in Eq. ([T|)] all the systems evolve independently 
according to their own dynamics. On increasing the cou- 
pling strength, the system xi{t) starts to drive the sys- 
tem X2{t). Consequentially, the system X3{t) starts to 
follow the drive system xi{t) for larger values of e and 
this is continued upto the N*'^ system. Hence, global 
synchronization is achieved via sequential synchroniza- 
tion of the systems in the array as a function of coupling 
strength. To be more specific, upon increasing the cou- 
pling strength, e, from zero, nearby systems to the drive 
in the array synchronize sequentially with it, while the 
faraway systems are still in their transition state. The 
other desynchronized systems will synchronize sequen- 
tially for further larger values of e. The occurrence of 
sequential phase synchronization in an array of unidirec- 
tionally coupled time-delay systems has been shown in 
Ref. [26| and sequential desynchronization in a network 
of spiking neurons is reported in Ref. [13] as a function of 
coupling strength, e. We may also note here a somewhat 
analogous situation occurs but now as a function of time 
for a fixed coupling strength in an array of unidirection- 
ally coupled chaotically evolving systems in Refs. (35l[36j . 

Here, we find that locally the synchronization in the 
array occurs immediately in a very narrow range of cou- 
pling strength; globally, it occurs in a broader range 
of e due to sequential synchronization. To character- 
ize these local and global synchronization transitions, we 
have calculated the probability of synchronization $i(£) 
(which is defined as the fraction of time during which 
|a;i(t) — a;i(i)| < S occurs, where (5 is a small but an arbi- 
trary threshold) and the average probability of synchro- 
nization [^'(e) = Si=2 ^«(^)]- Here the asynchro- 
nized state is characterized by $i(£) — 0, CS by $i(e) = 1 
and the transition region by intermediate values less than 
unity. 

To understand the dynamical organization of sequen- 
tial synchronization in the array (Eq. ([IJ), we have calcu- 
lated the probability of synchronization as a function of 
£ and the system index i, which is depicted in Fig. [1] 
In this figure, the black color indicates the asynchro- 
nized state {^i{s) = 0.0) and the yellow color (light 
gray) corresponds to the complete synchronization state 
(^^(e) = 1.0), while intermediate colors represent the 
transition region. From this figure one can clearly see the 
occurrence of sequential synchronization as a function of 
s where the nearby systems to the drive get synchronized 
first for lower values of e, whereas the far away systems 
are synchronized at larger e. 

Sequential synchronization can also be visualized using 
snap shots of the oscillators in the node vs node plots. 
We regard the oscillators in the array as synchronized 
when the probability of synchronization ^i{e) > 0.96, 
which are indicated by filled circles. Figure [5] shows 
node vs node diagrams for various values of the coupling 
strength. For e = 0.4 none of the oscillators are syn- 
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FIG. 4: The difference between some selected piecewise linear time-delay systems (Eq. ([T])) shows intermittent synchronization, 
(a) Aa;i.i5 for e = 0.76, (b) A2;i,2o for e = 0.81, (c) Axi.25 for e = 0.86, (d) Ax5,20 for e = 0.76, (e) Aa;5,25 for e = 0.81 and 
(f) Aa;5,30 for e = 0.86. 




FIG. 5: (Color online) The statistical distribution of the lam- 
inar phase for the systems i = 15 for e = 0.76, i = 20 for 
£ = 0.81 and the system i = 30 for e = 0.93, all satisfying — | 
power law scaling. 



FIG. 6: The average difference (AX) of all (A'^ — 1) piecewise 
linear time-delay systems in the array (Eq. Q) with the drive 
x\ shows an intermittent synchronization transition, (a) e — 
0.85, and (b) e = 0.89. 



chronized with the drive system (see Fig. Efa)). Figure 
IDJb) indicates that the first seven oscillators are synchro- 
nized with the drive for e = 0.7. Further increase in 
the coupling strength results in increase in the size of the 
synchronized cluster resulting in the formation of sequen- 
tial synchronization. Figures ^c) and [D^d) are depicted 
for e = 0.87 and 1.1, respectively, illustrating sequential 
synchronization. 

To discuss the nature of the synchronization transition 
locally, we have calculated the probability of synchro- 
nization for some selected systems (« — 15, 20 and 30) 
in the array as a function of e (see Fig. ^i{e) of the 
system i = 15 is plotted as a function of e (represented 
by the filled squares). In the range of e G (0,0.76), there 
is an absence of any entrainment between the systems re- 
sulting in an asynchronous behavior and $15 (e) is prac- 
tically zero in this region. However, starting from the 
value £ = 0.76 and above, there appear some finite val- 
ues less than unity attributing to the transition regime. 
Beyond e = 0.78, $15(5) attains unit value indicating CS. 
We have also plotted the probability of synchronization 
in Fig. 13] for two more selected systems i = 20 and 30 



represented by the asterisk symbol and filled triangles, 
respectively, indicating the immediate transition to CS 
locally. The system i = 20 attains CS at e = 0.84 and 
the system i = 30 reaches the CS state at e = 0.98. From 
this figure, one can understand the occurrence of sequen- 
tial synchronization of the individual systems (locally) in 
the array as a function of the coupling strength. To ex- 
plain global (macroscopic) synchronization phenomenon, 
we have calculated the average probability of synchro- 
nization (^'(e)) of the iV = 30 systems as a function of 
£ and depicted it in Fig. [3] (represented by the filled cir- 
cles). It confirms sequential synchronization by gradual 
increase in '3/(e) as a function of e (which indeed exactly 
matches with Fig.[T]). 

Next, in the transition regime, we observe intermittent 
synchronization in every individual system and this can 
be characterized qualitatively by a difference in the mag- 
nitudes of the states between the systems {Axi^i = \xi — 
Xi\,) for selected ones {i = 15,20 and 25). Figure UJ^a- 
c) shows intermittent synchronization in the above men- 
tioned systems for e = 0.76, 0.81 and 0.86, respectively. 
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FIG. 7: (Color online) The probability of synchronization 
"I>i(e) as a function of e and the system index i illustrating 
the occurrence of instantaneous synchronization transition to 
CS both locally and globally in mutually coupled piecewise 
linear time-delay systems (Eq. Q). 

We also find that the synchronization quality in the tran- 
sition region depends on the respective positions of the 
response systems from the drive, as well as on the dis- 
tance between the two units in the system and the cou- 
pling strength. We have additionally plotted the differ- 
ence between the systems A5^207 ^5.25 and As^so for the 
above set of values of coupling strength (Fig. IH^d-f)) to 
demonstrate the above features. 

The statistical features associated with the intermit- 
tent dynamics is also analyzed by the distribution of the 
laminar phases K{t) with amplitudes less than a threshold 
value of A (here we have choosen A = |a;i(t) — Xi{t)\ = 
1 0.001 1 , i — 15,20 and 25). A universal asymptotic 
power law distribution A(i) oc is observed for the 
above threshold value of A with the exponent a = —1.5. 
Fig. [5] shows the laminar phase distribution of the above 
selected systems. The filled circles represent a laminar 
distribution of the system z = 15 for e = 0.76, the filled 
triangles correspond to the laminar distribution of the 
system ? = 20 for e = 0.81 and the filled squares rep- 
resent a laminar distribution of the system z = 30 for 
e = 0.93 which clearly display the — | power law scal- 
ing, a typical characterization of on-off intermittency. It 
should be noted that this result does not change for a 
large range of A. 

To understand the phenomenon of intermittent syn- 
chronization transition globally in the whole array, 
we have calculated the average difference (AX — 
Y^^=2 1^1 ~ ^jl) ~ 1) systems in the ar- 

ray with the drive xi. Figures [6ja) and |6l[b) show the 
average difference for the coupling strengths e = 0.85 and 
0.89, respectively. 

The reason behind sequential synchronization transi- 
tion is in accordance with the sequential stabilization of 
all the unstable periodic orbits of the response systems 



FIG. 8: (Color online) Snap shots of node vs node plots of 
mutually coupled piecewise linear systems indicating instan- 
taneous synchronization, (a) e = 58.0, (b) e = 61.0. 

in the array as a function of the coupling strength. It 
is a well established fact that a chaotic/hyperchaotic at- 
tractor contains an infinite number of UPOs of all pe- 
riods. Synchronization between the coupled systems is 
said to be stable, if all the UPOs of the response sys- 
tems are stabilized in the transverse direction to the syn- 
chronization manifold. Consequently, all the trajectories 
transverse to the synchronization manifold converge to 
it for suitable values of e. For sequential synchroniza- 
tion, the UPOs in the complex synchronization mani- 
fold of the response systems near to the drive are sta- 
bilized first for appropriate threshold values of the cou- 
pling strength e as it is increased, while the UPOs of the 
far away systems remain unstable for these values of e. 
Once the coupling is increased further, the UPOs of the 
far away systems are gradually stabilized as a function of 
the coupling strength. Unfortunately, methods for locat- 
ing UPOs have not been well established for time-delay 
systems, which has hampered a qualitative proof for the 
gradual stabilization of UPOs by locating them. 

III. SYNCHRONIZATION IN A PIECEWISE 
LINEAR TIME-DELAY SYSTEMS: LINEAR 
ARRAY WITH BIDIRECTIONAL COUPLING 

In this section, we consider an array of mutually cou- 
pled (bidirectional coupling) piecewise linear time-delay 
systems with identical subunits. The dynamical equation 
then becomes 

Xi = -axi{t) + /3 f {x,,{t-T))+e[xt+i{t) - 2xi{t) +Xi^i{t)], 

(4) 

where i = 1,2, ••• ,iV. We choose open end bound- 
ary conditions. The parameter values are the same as 
in Sec. [Ill The nonlinear function f{x) is chosen as in 
Eqs. (12) -(O. In the mutual coupling case there are no 
drive and/or response systems where each and every os- 
cillator shares the signals mutually with its two nearest 
neighbors. So the synchronization transition is instan- 
taneous due to the mutual sharing of the signals and 
one needs a very large value of e to attain CS. In the 
transition regime, we have observed an intermittent syn- 
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FIG. 9: (Color online) The probability of synchronization 
$i(e) of selected systems {i = 15,20 and 30) and the av- 
erage probability of synchronization ('l'(e)) in the mutually 
coupled array (Eq. (|4])) as a function of e. 

chronization transition in all the systems simultaneously 
in the array. 

We have calculated the probability of synchronization 
of all the = 30 systems in the array as a function of the 
coupling strength e and the system index i (see Fig.[71). In 
this figure the black color represents the desynchronizcd 
state ($(e) = 0.0) and CS is represented by the yellow 
color (light gray) ($(e) = 1.0). The transition regime is 
indicated by intermediate colors. From this figure one 
can clearly see that locally every individual system re- 
quires large values of e to attain CS and globally all the 
systems synchronize immediately for the same value of e. 

In the mutually coupled array, all the systems get syn- 
chronized immediately in a narrow range of the coupling 
strength, in contrast to sequential synchronization. Fig- 
urelD^a) is plotted for e = 58.0, where none of the oscilla- 
tors in the array are synchronized, whereas for e = 61.0, 
all the systems are completely synchronized as depicted 
in Fig. lib). 

To characterize the nature of synchronization transi- 
tions to CS both locally and globally, we again use the 
probability of synchronization $(£) and the average prob- 
ability of synchronization 4'i(e), respectively. In Fig. [HI 
we have plotted <&,(£) for some selected piecewise lin- 
ear systems (i — 15, 20, 30) as a function of e. For in- 
stance, we have illustrated ^i{e) for the system i = 15 
in Fig. [9] (represented by the filled squares). From this 
figure, one can observe that in the range of e G (0, 50) 
there is an absence of any entertainment between the 
systems resulting in asynchronous behavior and $15 (e) 
is low ($i5(e) < 0.4). However, for e > 50 there appear 
oscillations in $15 (e) in the range of e e (50, 60) exhibit- 
ing intermittent transition. Beyond e = 60.0, $15 (e) = 1 
indicating perfect CS of the system i — 15. We have 
also calculated $(e) for the systems i = 20 and i = 30, 
represented by asterisk symbols and filled triangles, re- 
spectively, which show similar transitions to CS almost 
at the same value of e. We have also confirmed a simi- 
lar immediate transition to CS in all the systems in the 
array (see Fig. [7|) . 

To examplify the global synchronization phenomenon. 



we have calculated the average probability of synchro- 
nization (^'(e)) of TV = 30 systems as a function of the 
coupling strength as shown in Fig. |9]by the filled circles. 
In the range of e S (0,54), there is an absence of any 
synchronization and so ^'(e) is having zero or low values 
(^'(e) < 0.3). In the range of e e (54, 59), ^'(e) is charac- 
terized by some finite values less than unity and beyond 
e > 59.0 there is a sudden jump to the value of ^(e) = 1.0 
corrobrating all the systems synchronize immediately at 
the same value of the coupling strength attributing to 
the occurrence of global CS. Further, in the transition 
region we have found long time intervals of high quality 
synchronization which is interrupted at irregular time in- 
tervals by intermittent chaotic bursts simultaneously in 
all the systems in the array which we call as CIS. 

To demonstrate the existence of CIS, we have calcu- 
lated the spatiotemporal difference {Ax{t) = \xi{t) — 
Xi{t)\, i =^ 2,3, ••• ,30) of the array as a function of 
time and the oscillator index i as in Fig. [TU] for different 
values of e. Here the black color indicates zero differ- 
ence {Ax{t) — 0.0) and the red and green colors (dark 
and light gray) indicate bursting amplitudes. In the ab- 
sence of coupling (e = 0.0), the systems are evolving 
independently and so there is no correlation between the 
systems as shown in Fig. fTOFa) . If we increase the cou- 
pling, we observe several intermittent bursts along with 
the synchronization as depicted in Figs. iroT b) and iroT c) 
for e = 50.2, and 54.6, respectively. From these figures, 
one can clearly see the occurrence of aperiodic intermit- 
tent chaotic bursts along with the synchronized regions 
simultaneously in all the systems in the array. Beyond 
e > 60 one can observe CS as illustrated in Fig. [TUT d) 
where the spatiotemporal difference of the systems is ex- 
actly zero for e = 61.0. 

To elaborate the occurrence of CIS in the array more 
clearly, we have calculated the difference between the 
systems coupled in the array and plotted for some se- 
lected systems {i = 20 and 30) for different values of 
the coupling strength. The difference between the sys- 
tems 1 and 20, Axi^2o{t) = \xi{t) — X2o{t)\, is plotted for 
e = 50.2, 54.6 and 58.1 in Figs.IIIta) -^c), respectively, 
which clearly displays the existence of aperiodic intermit- 
tent bursts along with the synchronized regions. We have 
also plotted the difference Aa;i,3o(t) = |a;i(i) — a:3o(t)| for 
the same values of e as shown in Figs. fTlT e) - fTlT g). It 
is to be noted that in both systems {i = 20, 30) the in- 
termittent bursts simultaneously occur at the same time 
and this occurs in all the other systems connected in the 
array confirming the existence of CIS. In Figs. ITlT d) and 
[TlT h) the difference between the systems completely van- 
ish for £ = 61.0 indicating the occurrence of CS. We have 
plotted the above figures with 10''' time units after leaving 
a sufficient number of transients. 

Further statistical features associated with the inter- 
mittent dynamics of the entire array are also analyzed by 
calculating the distribution of the laminar phases A{t), 
which is shown in Figs. Il^ a) and[T2jb) for selected sys- 
tems i = 20 and 30, respectively, for e = 54.6 and 58.1 



FIG. 10: (Color online) The spatiotemporal difference {Ax(t)) of the mutually coupled piecewise linear systems for various 
values of coupling strengths, (a) e = 0.0, (b) e — 50.2, (c) e — 54.6, (d) e = 61.0. Here the black color indicates that the 
difference is zero and the red and green colors (dark and light gray) indicate the bursting amplitudes. 
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FIG. 11: (a)-(d) The difference between systems 1 and 20 (Aa;i,2o = \xi{t) - X2o{t)\) for e = 50.2, 54.6, 58.1 and 61.0. (e)-(h) 
The difference between systems 1 and 30 (Axi.ao = \xi{t) — X3o{t)\) is plotted for the same set of coupling strength values 
given above. 




rectional coupling configuration is due to the fact that 
the strength of the couphng e contributes only less sig- 
nificantly to stabilize the UPOs as the error in the cou- 
pling term in Eq. Q gradually becomes smaller from the 
transition regime after a certain threshold value of the 
coupling strength. 



FIG. 12: (color online) The statistical distribution of the 
laminar phase for selected piecewise linear systems (Eq. Q 
(a)j = 20 and (b) i = 30 satisfying a — | power law scaling 
for the coupling strength e = 54.6 and 58.1. 



which clearly display the — | power law scaling to confirm 
the on-off intermittency. 

The reason for the occurrence of GIS can be explained 
as follows: As we have already explained, a chaotic at- 
tractor can be considered as a pool of infinitely many 
UPOs of all periods. Synchronization between the sys- 
tems are asymptotically stable, if all the UPOs of the 
systems are stabilized in the transverse direction to the 
synchronization manifold. Consequently, all the trajecto- 
ries transverse to the synchronization manifold converge 
to it for suitable values of the coupling strength and this 
is reflected in the stabilization of the UPOs upon syn- 
chronization. From our results, we find that the UPOs of 
the systems are stabilized in the complex synchronization 
manifold only for a very large value of coupling strength 
after a certain threshold value. It is also to be noted 
that the intermittency transition in the case of a bidi- 



IV. CONCLUSION 

In conclusion, we have shown the existence of sequen- 
tial and instantaneous synchronization transitions in an 
array of time-delay systems with different coupling con- 
figurations. If the systems are coupled with unidirec- 
tional configuration, we have observed an immediate syn- 
chronization transition to CS microscopically and if we 
consider the macroscopic synchronization behavior of the 
entire array we find that the transition region is gradu- 
ally increasing as a function of e due to sequential syn- 
chronization which is verified by the probability of syn- 
chronization and average probability of synchronization. 
In the transition regime we have observed the existence 
of intermittent synchronization. On the other hand, if 
we consider an array of mutually coupled time-delay sys- 
tems, every individual system (microscopically) synchro- 
nizes immediately for a very large value of s and globally 
(macroscopically) the synchronization transition occurs 
immediately in the whole array. In the transition region 
a new type of synchronization called GIS occurs which is 
characterized by long intervals of high quality synchro- 
nization interrupted at irregular times by intermittent 
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chaotic bursts simultaneously in all the systems. 

The reason (mechanism) for these two distinct tran- 
sition scenarios is explained based on unstable periodic 
orbit theory. The GIS is confirmed using the spatiotem- 
poral difference and a power law behavior of the lami- 
nar length distributions with — | power law scaling. The 
above studies have been carried out in a well known 
piecewise linear time-delay system. We have also con- 
firmed the occurrence of the above results for another 
well known time-delay systems namely the Mackey-Glass 
system [s^ with an array length of = 50 and we do 
observe the same kind of sequential and instantaneous 
synchronization transitions preceded by GIS for unidirec- 
tional and bidirectional coupling configurations, respec- 



tively. 
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